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ABSTRACT 

We present age measurements for 195 star clusters in the Small Magellanic Cloud based on com- 
parison of integrated colors measured from the Magellanic Clouds Photometric Survey with models 
of simple stellar populations. We find that the modeled nonuniform changes of cluster colors with age 
can lead to spurious age peaks in the cluster age distribution, that the observed numbers of clusters 
with age, t, declines smoothly as t~ 2A , that for an assumed initial cluster mass function scaling as 
M ~ 2 the dependence of the cluster disruption time on mass is proportional to M 0AS , that despite the 
apparent abundance of young clusters the dominant epoch of cluster formation was the initial one, and 
that there are significant differences in the spatial distribution of clusters of different ages. Because of 
limited precision in our age measurements, we cannot address the question of detailed correspondence 
between the cluster age function and the field star formation history. However, this sample provides 
an initial guide for which clusters to target in more detailed studies of specific age intervals. 
Subject headings: globular clusters: cluster ages — galaxies: evolution — galaxies: individual (Small 
Magellanic Cloud) — galaxies: star clusters — Magellanic Clouds 



1. INTRODUCTION 

The dis covery of young star cl usters in ongoing galaxy 
mergers ijHoltzman et all Il992t iWhitmore et all 119931) 

has provided added motivation for understanding the 
formation mechanism of stellar clusters and the dynam- 
ical processes that lead to their destruction. While the 
ancient clusters in our own galaxy provide only indi- 
rect information regarding the formation or disruption 
of clusters, the Magellanic Clouds cont ain numerous 
young clusters (< 1 Gyr old) for study ijHodgd 119611: 
Ivan den"Be rgh 198.3). The large number of clusters is of 
great statistical value, but it has divided the studies into 
two camps. Ages have been determined either for large 
samples of clusters using low-precision measurements , 
such a s integrated colors (for example see Ivan den B orah 
l|1981l) ) or for small samples using high- precision mea- 
surements, such as color-magnitude diagrams (for an 
early example see iBaird et all ()1974|) 'l. Each approach 
has its strengths and weaknesses, but a study of the 
overall population must include age estimates for a large 
fraction of all the clusters in the galaxy. Only with such 
measurements can one begin to address whether the clus- 
ter formation history tracks the field star formation his- 
tory, how quickly clusters are dissolved, and whether the 
cluster system as a whole still retains some memory of 
cluster formation episodes. 

The study of the clusters in the Magellanic Clouds 
has a rich history that cannot be properly summarized 
here. Nevertheless, from that work we know of several 
interesting features in the age distribution of clusters in 
both the Large and Small Magellanic Clouds (LMC and 
SMC). In the LMC there is an "age gap" (a deficit of clus- 
ters of ages corresponding to ages between roughly 3 and 
10 Gy r) that ha s been extensively studied llensen et al " 

)]j IRich 



et al 



19881 IDa Costal Il99lt Ivan den Berehl IT99 

2001). The age distribution of clusters in the SMC 
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less established. Ground based studies of SMC clus- 
ters hav e concluded that the age distribution is con- 
tinuo us l)Da Costa fc Hatzidimitrioulll998tlMighell et alJ 
1998), but analysis of the SMC's seven brightest old (age 
> 1 Gyr) clusters with the Hubble Space Telescope has 
suggested two distinct episodes of clus ter formation tha t 
occurred 2 ± 0.5 Gyr and 8±2 Gyr ago <|Rich et alte OOOL 
Three additional SMC clusters that satisfy the age crite- 
ria but are too faint to be included in the HST study are 
Lindsay 1, which has a ground-based age of 9 Gyr and so 
is coincident with the older burst, Lindsay 11, which has 
an age o f 3.7 Gyr, and Lindsay 113, which has an age of 
5.3 Gyr ijMighell et allll99^ IRich et a1.ll27jonJl Neither 
of the latter two appear to have formed within the two 
identified bursts and may populate the "gap" . Whether 
the structure in the age distribution is real, and whether 
it tracks any pattern observed in the field star formation 
history of the SMC, are open questions. 

We present a study based on integrated colors of 195 
clusters in the SMC that incorporates a number of im- 
provements over similar previous studies: 1) we use dig- 
ital imaging in four filters (most previous studies of in- 
tegrated colors were based on photographic plates and 
utilized only two or three filters), 2) we use the struc- 
tural parameters in dividually derived for each cluster 
( Hill & Zaritsk vl2003|) to set our aperture size (most pre- 
vious studies have a fixed aperture size for all clusters), 
and 3) we use multiple theoretical models of simple stel- 
lar populations to derive an age and associated uncer- 
tainty (most previous studies calibrated integrated col- 
ors to a small number of clusters with ages derived from 
color-magnitude diagrams — a procedure that can lead 
to coarse age resolution and unknown systematic prob- 
lems that are difficult to explore, see §3.1). The data and 
methodology used to measure ages are described in §2. 
We discuss various aspects of the derived age distribution 
in §3 and summarize in §4. 
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Fig. 1. — UBVI images of the clusters. Each row contains the 
four filter images for two clusters. The images are 70 arcsec on a 
side and oriented such that N is to the right and E to the bottom. 
This Figure includes the images for clusters 1 through 16. The 
remainder of the images are available in the electronic version. 



2. DATA AND ANALYSIS 
2.1. Images of the Cluster Sample 

The sa mple of SMC clusters that we analyze is pre- 
sented bv lHill fe Zaritskv! l)2003j) . Coordinates and cross- 
identifications are given in Table 2] With the exception 
of clusters and associations embedded in the complex en- 
vironment of star forming regions, these 204 clusters are 
the most unambiguous clusters in the 4.5° x 4° section 
of the SMC o bserved by the Mage llanic Clouds Photo- 
metric Survey l|Zaritskv et alJl2002j) . Figure ^ shows im- 
ages of clusters 1 through 16. We have excluded clusters 
and associations that lie within emission line nebulae be- 
cause of the difficulty in measuring integrated colors in 
such an environment. The absence of these clusters ex- 
plains in part why other catalogs list man y more star 
clusters (for exampl e see iPietrzvhski et all l|1998fl and 
iBica fe Dutral pOOO^ . but wc also found from visual in- 
spection of over a thousand candidate clusters that many 
are quite marginal. It is not our contention that only 
the clusters presented here are real. There are almost 
certainly real clusters lurking among the previously pub- 
lished catalogs that we have excluded. However, even if 
real, we cannot confidently measure the colors or struc- 
tural parameters of such poor clusters. When consider- 
ing the results from our sample, one must remain aware 
that we are incomplete in both the youngest and poorest 
(poor in either total number of stars or surface density) 
clusters. The latter might be a particularly interesting 
population because it may consist of clusters that are 
nearing the end of their lifetime as bound objects. 

Our cluster images are drawn from drift scan imaging 
of the central 4° x 4.5° of the SMC in U, B, V, and / 
done with the Las Campa nas Swope telescope (lm) and 
the Great Circle Camera l|Zaritskv et al.lll996f) between 
1996 November and 1999 December. The effective expo- 
sure time for any portion of the SMC is between 4 and 5 
min and the pixel scale is 0.7 arcse c pixel -1 . A photomet - 
ric catalog of stars is presented by Zaritsk v et al.l (2002), 
but in this study we use the reduced images rather than 
the stellar catalog because of the catalog's incomplete- 
ness in the high-density centers of stellar clusters. We 



extract 350" x 350" subimages in each of the four fil- 
ters centered on the cluster from the larger drift scans. 
We use the photometric solutions that were applied to 
the stellar catalog, which have an observational scatter 
of between 0.01 and 0.04 mag for standard star fields. 
This uncertainty is significantly smaller than that result- 
ing from other sources of error in the measurement of 
cluster integrated colors. 

2.2. Integrated Cluster Colors 

The measurement of an integrated color or magnitude 
depends critically on the degree to which one can sub- 
tract the effect of the underlying background. This is 
one area in which the availability of moderate resolution 
digital images should enable a quantitative improvement 
over previous studies. The two primary considerations 
are where to set the cluster and background apertures 
and how to remove contamination by stars that are too 
bright to be plausible cluster members or random SMC 
field stars (these are typically foreground Galactic stars). 

The optimal choice of aperture size is not evident. The 
clusters in this sample vary significantly in size, so a sin- 
gle aperture size is inappropriate. Although the cluster's 
tidal radius might be a natural choice for the outer scale 
of the aperture, t he measurement of th e tidal radius is 
highly uncertain. Hill & Zaritskv (2003) found that the 
radius that encloses 90% of the light of a cluster, r go , is 
much less sensitive to uncertainties in defining the back- 
ground level. We use their tabulated values of r 90 in the 
V band to set our aperture size for each cluster indepen- 
dently. The apertures are circular and are the same for 
all filter bands. The background is calculated from the 
area of the 350" x 350" cluster image that lies 5 pixels 
(3.5") beyond r 90 . 

The effect of contamination in both the cluster and 
background apertures can be severe. We minimize the 
impact of the brightest stars on our calculation of the 
background by adopting an upper limit to the pixel val- 
ues used in the calculation. First, we calculate the mean 
background, without this threshold. Then we set the 
threshold to be equal to the mean background level plus 
two times the central flux value of the cluster. Lastly, we 
calculate the mean background level using this threshold. 
The possibility of foreground contamination within the 
cluster aperture is smaller because its area is smaller than 
that of the background aperture. However, from visual 
examination we suspect that the photometry of as many 
as 73 clusters may be significantly contaminated by an 
unusually bright star for at least one of the filters within 
the aperture corresponding to the upper la bounds on 
rgo. Because the measurement of the cluster luminosity 
is the integrated flux within the aperture, rather than 
an average, we cannot simply exclude high valued pixels. 
Instead we flag such clusters in Tabled but do not cor- 
rect their photometry. We found that these clusters do 
not have noticeably different properties than the others. 

To arrive at the measurement of the cluster magnitude 
in each filter, we sum the counts within the rgo cluster 
aperture, subtract the mean background level, multiply 
by 10/9 to correct for the choice of rgo as the aperture, 
and iteratively apply the photometric solution to incor- 
porate the color terms. Magnitudes and colors are pre- 
sented in Table ^ To estimate uncertainties in these 
quantities, we recalculate the magnitudes and colors us- 
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TABLE 1 
CLUSTERS AND THEIR COLORS 



Number RA DEC M v a v U-B cr u _ B B-V cr B _y V-I <j V -i Name" 

1 24 18.67 -73 59 35.8 15.92 1.17 0.75 0.24 0.72 0.14 0.92 0.07 

2 24 43.16 -73 45 11.7 K5, L7, ES028-18 

3 25 26.60 -74 4 29.7 14.23 ■■■ 0.57 ■■■ 0.65 ••■ 1.10 ••■ K6, L9, ESO28-20 

4 27 45.17 -72 46 52.5 13.95 0.38 0.03 0.09 0.80 0.04 1.20 0.05 K7, Lll, ES028-22 

5 28 2.14 -73 18 13.6 14.85 1.17 0.28 0.05 0.61 0.05 1.00 0.40 K8, L12 

6 29 55.22 -73 41 57.1 14.31 1.12 -0.05 0.48 0.65 0.06 0.37 0.32 HW3 

7 30 0.26 -73 22 40.7 13.84 ■■• 0.54 ■■• 1.08 ■■■ 1.32 ■■• K9, L13 

8 31 1.34 -72 20 30.0 14.64 0.28 0.00 0.04 0.80 0.07 1.20 0.12 HW5 

9 32 41.02 -72 34 50.1 14.66 0.35 0.26 0.18 0.87 0.05 0.90 0.05 L14 



Note. — The complete version of this table is at: http://ngala.as.arizona.edu/dcnnis/rafclski.html The printed edition contains 
only a sample. 

a K is for lKronl 1195611. L is for | Lindsavl 119581). HW is fo r lHodge fc Wrighti tl97l . H is for lHodgel 1198511 . H86 is for lHodgel I1986D . 
BS95 is for lBica fc SchmitH I1995h . B is for lBriicll 1197611 



ing the la bounds on rgr,. iHill fc Zaritskvl l)2003|) do not 
present oy go for some clusters because their best fit King 
model was not statistically acceptable. For these clusters 
we adopt the average of the color uncertainty calculated 
for the other clusters. The uncertainty from the pho- 
tometric calibration is added in quadrature, but it is a 
minor contribution. Magnitudes and colors, in particu- 
lar filter bands, are unavailable for the 35 clusters whose 
images are saturated in at least one of the filters, the 
two clusters that are off an edge of a scan, and the one 
cluster (number 180) whose image quality is too low to 
be useful. These problems often affect the cluster image 
in only one filter, and in such cases the magnitudes are 
calculated for the unaffected filters. There are 38 clusters 
that are missing at least one of the three colors in our Ta- 
ble ^ but only seven for which we are unable to compute 
any colors. Eliminating these seven plus the two that are 
off the edges of scans leaves us with 195 clusters, out of 
the original 204, for which we measure an age. 

Because contamination is such a serious problem, and 
because, as we will show below, the integrated colors 
scatter widely around the model predictions, we investi- 
gate whether adopting a smaller cluster aperture would 
lead to more robust colors. Using colors from small aper- 
tures assumes that there are no internal color gradients. 
We re-measure colors using an aperture of radius 0.5rgo 
and find no noticeable decline in the scatter of the colors 
about the model predictions. We conclude that the con- 
tamination cannot be significantly mitigated by decreas- 
ing the aperture, partially because the scatter is also due 
to stochastic effects within the cluster's own population. 

We obtain an external estimate of the measurement 
uncertainties by comparing our magnitudes and colors 
to previous studies where possible. A compilation of in- 
tegrated c olors from various ph otometric studies is pre- 
sented by Ivan den Berghl itlQSl . Figures [5] and |3 show 
the differences bet ween our V magnitude s and colors and 
those tabulated bv lvan den Berghl l)1981|) . The error bars 
are underestimated because they do not include the un- 
certainties in the compiled data (none were quoted) . As- 
suming a standard error of 0.1 mag, the scatter in B — V 
appears to be entirely consistent within the uncertain- 
ties, while it appears that we have slightly underesti- 
mated the uncertainties in U — B (either in our data or 
van den Bergh's). In the mean both colors agree well. 
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Fig. 2. — Co mpariso n of our V m agnitudes, VrZ^ with those from 
Ivan den Berghl <198ll) . V v dB. The lvan den Berghl ll9olT) compila- 
tion did not quote photometric uncertainties, and thus only uncer- 
tainties in our measurements are shown. The errorbars correspond 
to those in the tables. The line is the 1:1 correlation. 



To quantitatively determine whether the scatter is com- 
mensurate with the quoted uncertainties, we examine the 
distributions of Acolor/er. If the uncertainty estimates 
are correct, a Gaussian fit to the distribution should have 
cr=l. Instead, we find (after removing the mean off- 
set, one highly discrepant cluster, and not adopting any 
uncertainties for the published data) that the fit has a 
dispersion that is ~ 1.5. However, if the uncertainties 
in the published data are similar to ours, then cta would 
be 40% larger than calculated and the fitted Gaussian 
would have a dispersion that is close to one. We con- 
clude that except for the offset in the mean, our colors 
are statistically consistent with previous measurements 
and our uncertainty estimates are appropriate. 

2.3. Determining Cluster Ages 

The integrated colors of any stellar population provide 
a luminosity weighted age measurement. The color 
should be a particularly stable chronometer for clusters 
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Fig. 3. — Comparison of our B — V and U — B colors with 
those from van den Berghl 119811). The differenc e is our values - 
Ivan den Berghl 119811) . The Ivan d en Bcrgh (1981) compilation did 
not quote photometric uncertainties, and thus only uncertainties in 
our measurements are shown. The errorbars correspond to those 
in the tables. 



because they formed their stars over a timescale that is 
much shorter than their age. We can measure the age of 
a cluster to the degree that we can accurately measure 
and model the colors. The colors are pre sented in Table 
|T1 and the models we use are those of iLeitherer et al.l 
(1999), hereafter ref erred to as the Starburst99 or 
S99 m odels, and of I Anders and Fritze-v. Alvenslebenl 
(2003), hereafter referred to as the GALEV mod- 
els. The S99 models that we use correspond to the 
standard mass loss prescription, the theoretical wind 
model, and the full isochrone mass interpolation. We 
have adopted a Salpeter IMF in both models. We 
derive ages using models with a range of appropriate 
metallicites (Z = 0.001, 0.004, and 0.008 for the S99 
models and 0.004 and 0.008 for the GALEV models). 
These metallicity values refer to the mass fraction of 
metals relative to hydrogen and correspond to [Fe/H] 
of —1.3, —0.7, and —0.4, respectively. Observations of 
individual clusters (iDopita et all 119851 iDa Costal 119911: 
Page! fc Tautvaisienei |199*8[ Ide Freitas Pacheco et al.l 
19981 iPiatti et al.l 120011) and the reconstruction of the 
global chemical enrichment history of the SM C from 
color-magnitude diagrams l)Harris fc Zaritskvl |2004() 
suggest that the metallicity of the SMC has for the 
most time varied between [Fe/H] ~ —1.2 and —0.4. 
As we discuss below, some of the finer details of the 
cluster formation history are extremely sensitive to 
the exact choice of metallicity for the model and our 
current constraints on the age-metallicity relationship 
are insufficiently precise to dictate which model should 
be used for each age. 

In Figures^ andEJwe overlay the two sets of models 
over our reddening-uncorrected cluster colors. We find a 
mean agreement with the S99 models, although the ob- 
servational scatter is large. The principal systematic dif- 
ference between our data and the models appears to be 
in the V — I colors (which are too red in the models for 
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Fig. 4. — Integrated U — B vs B — V from data and models 
for SMC clusters. The triangles represent our measurements, the 
points, which blend into a line in certain regions, represent the 
Starburst99 model for Z = 0.004 (for ages between 0.001 and 10 
Gyrs), and the line represents the GALEV model for Z = 0.004 
(for ages between 0.004 and 14 Gyrs). 
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Fig. 5. — Integrated V — I vs B — V from data and models 
for SMC clusters. The triangles represent our measurements, the 
points, which blend into a line in certain regions, represent the 
Starburst99 model for Z = 0.004 (for ages between 0.001 and 10 
Gyrs), and the line represents the GALEV model for Z = 0.004 
(for ages between 0.004 and 14). 



the bluest clusters). The general agreement of the B — V 
colors suggests that the problem lies in the modeling of 
the /-band. The agreement is significantly poorer rela- 
tive to the GALEV models, but we continue to use those 
models as well to examine systematic errors introduced 
into the derived ages by the use of one or the other set 
of theoretical models. 

The observed scatter of the colors about the theoreti- 
cal model is the fundamental observational limit to the 
precision with which we can measure ages. The scatter 
arises from a variety of sources including the lack of a 
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Fig. 6. — Integrated U — B vs V — I from data and models 
for SMC clusters. The triangles represent our measurements, the 
points, which blend into a line in certain regions, represent the 
Starburst99 model for Z = 0.004 (for ages between 0.001 and 10 
Gyrs), and the line represents the GALEV model for Z = 0.004 
(for ages between 0.004 and 14). 



Fig. 7. — Comparison of ou r color measurements (triangles) to 
those of lvan den B crgh (1981) (circles). Straight lines connect the 
measurements for the same cluster. The models are as in Figures 
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reddening correction, which we discuss below, stochastic 
effects within the cluster population, contamination of 
the cluster apertures, which we discussed in §2.2, and 
the contamination of the background aperture, which 
we attempt to minimize. The last three of these pos- 
sibilities will become less important as we consider in- 
trinsically brighter clusters. A comparison (Figure 0) 
bet ween our U — B, B — V colors and those presented 
bv Ivan den Ber gh (1981), shows that indeed for these 
clusters, which are among the most luminous clusters, 
the scatter about the model line (in particular the S99 
model) is much smaller than for the sample as a whole 
(Figure 0}. We note, however, that the difference be- 
tween the two models is often comparable to the scatter 
in the data. 

To measure a cluster's age, we calculate the color dif- 
ference between data, typically U — B, B — V and V — I, 
although not all colors are available for all clusters (see 
§2.2), and models, at all ages available for the particular 
model. At each age, x 2 is evaluated using all the available 
colors and we adopt the age that corresponds to the mini- 
mum x 2 ■ Age limits corresponding to the 90% confidence 
range are derived using the criterion that A% 2 < 2.7. 
The best fit ages and uncertainties are given in Table 
with the subscript U for reddening-uncorrected and C 
for reddening corrected (see below for discussion of the 
reddening correction). 

Although the age fitting is a maximum-likelihood tech- 
nique, we also use the minimum xt to evaluate the good- 
ness of the best-fit model. There are a large number 
of clusters for which the model fit can be rejected with 
high confidence (for example, for the GALEV model with 
metallicity of 0.004 over half of the fits could be rejected). 
These clusters are not necessarily those with poor struc- 
tural fits or with obviously questionable photometry (for 
example, they are not those that appear to be contam- 
inated by a nearby bright star). Because the models 



do not account either for stochastic effects at the top 
end of the stellar luminosity function nor for contamina- 
tion, we suspect that these poor fits reflect such problems 
rather than some inability of the models to reproduce the 
cluster stellar populations. Supporting evidence for this 
conjecture comes from the similarity in the derived age 
distribution for clusters with low and high values of the 
best-fit xV s - We include all clusters, regardless of their 
minimum xt value, in subsequent discussion. 

One potential cause of high xt values is the omission, 
so far, of an extinction correction. We explore two op- 
tions for extinction corre ctions (both adopt a standard 
Galactic extinction law l)Schildl[i977[) which is accept- 
able for the SMC at optical wavelengths): 1) we in- 
clude extinction as a free parameters in our fitting algo- 
rithm, and 2) we adopt extinction values based on other 
data. Option 1 failed to produce reliable estimates of 
the extinction. Because of insufficient observational con- 
straints, the algorithm corrected for color scatter in a 
random way producing both brightening and dimming 
effects with larger corrections for more scattered clus- 
ters. We attribute this failure to that r ather small ex- 
tinction values expected across the SMC llZaritskvlll999D 
and to the degree of noise in our color measurements. In- 
stead, we choose to use th e extinction distributions mea- 
sured by IZaritskvl lj!999|) using thousands of individual 
stars across the SMC. The extinction is small along most 
lines of sight, but does vary as a funct ion of stellar type. 
We ad opt the approach presented bv lHarris fc Za ritskv 
l)2004[) where objects younger than 10 Myr are assigned 
the extinction derived from the young stars, objects older 
than 1 Gyr are assigned the extinctions derived from the 
older stars, and objects between these ages are assigned 
an extinction that is a linear i nterpolation (over l o g age) 
between the two boundaries. lHarris fc Zaritskvl (2004) 
assign objects a random extinction drawn from the ob- 
served distribution of extinctions at each location, here 
we simply assign the mean of the distribution. The ex- 
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TABLE 2 
AGES 



Number Agey Agejj — Ageu+ Agec Agec — Agec+ Contamination 



1 


1264 


1104 


2352 


1264 


976 


2420 




2 
3 


2520 


1368 


8080 


2148 


1160 


11500 




1 


6280 


5140 


7640 


5940 


3040 


7500 




5 


2080 


1232 


2920 


1264 


1136 


2220 


X 


6 


2288 


1264 


4000 


1996 


1148 


3760 




7 


12020 


7480 


14000 


12020 


5500 


14000 


X 


8 


924 


724 


1144 


316 


292 


532 




9 


4060 


2660 


6600 


2740 


1996 


6360 





NOTE. — Ages arc in Myrs. Agejj represents uncorrected Ages, and Agec present extinction 
corrected ages. An X under contamination denotes a contaminated cluster. The complete 
version of this table is at http://ngala.as.arizona.edu/dennis/rafclski.html along with other 
similar tables for different models and mctallicitics. The printed edition contains only a sample. 



tinctions are calculated and added to the model colors, 
for which the age, and hence corresponding extinction, is 
known. Comparison of x 2 values proceeds as previously 
to provide a best-fit measurement of the age. Although 
extinction should be included, it does not significantly 
decrease the scatter of colors about the models. 

We now compare our derived ages to those pre- 
sented in the literature. The results shown in Fig- 
ure |H1 include both a comparison to ages obtained 
via in tegrated colors <|van den Bergblll 981t iHunter et all 
2QQ3D and via iso chrone fitting l| Piefazvnslri et^lJ^^St. 
de Oliveira eTaTl l200?l iMierhell et all 119981: iRich et all 
2000). Ivan den Berghl l)1981[) provides only an age range 
for each cluster and we adopt the midpoi nt of the range 
for th e comparison. For the ages from Mighcll et al. 
(119981) we adopt the primary stated ages, and for the 
IRich et alJ IpOOOj) ages we used their [Fe/H] of -0.71 
ages that take into account the SMC distance modu- 
lus and assume color shifts are due to reddening. Our 
ages are broadly correlated to those presented in the lit- 
erature. However, when comparing to specific data sets 
differe nt patterns emerge. F or example, in comparison 
to the Ivan den Berghl l|1981|) data we appear to system- 
atically overpredict the ages. Ignoring any systematic 
differences, we calculate that the dispersion about the 
1:1 line for all the samples is 0.76. This result suggests 
that the ages are good to a factor of two. However, com- 
paring only to the most precise ages, those derived from 
color-magnitude diagrams, the scatter drops to 0.49, pro- 
viding added confidence in our measurements. 

3. DISCUSSION 

3.1. The Cluster Age Distribution for Ages < 5 Gyr 

We present histograms in Figure of the age distri- 
bution of stellar clusters in the SMC. Depending on the 
model, there are between 6 and 48 clusters with ages 
greater than 5 Gyr that are not shown in this Figure 
(but are discussed later). The Figure illustrates two key 
points: 1) independent of metallicity and model type, 
the current age distribution of clusters is strongly peaked 
to ages < 1 Gyr (this peak would be further enhanced 
by the inclusion of clusters embedded within emission 
line regions), 2) independent of metallicity and model 
type, the distributions are highly variable, although the 
location of peaks is strongly sensitive to metallicity and 
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Fig. 8. — Comparison of our measured cluster ages to those 
available in the literature with our values being Age^Z and the 
diff erent literature ages a re Age; it. Th e various literature studi es 
are Ivan den Bcrgn] 119811) 's (tr iangles). IPiet rzvnski et al. Il999)'s 
(filled bo xes). IRich et ^112000) (open boxes). IMighell et al l 119981) 
(cross es), IHunter et ail 120031) (filled circles) and lde Uliveira et alJ 
(2000)'s (open circles). For clarity, error bars are omitted. The 
line is the 1:1 correlation. 



model type. 

Although a histogram is one way of presenting the 
data, it fails to convey a sense of the age measurement 
uncertainties, which differ greatly among clusters. We 
present the smoothed age distribution in Figure 1101 in 
which each cluster is modeled by a normalized proba- 
bility distribution that is Gaussian on either side of the 
best-fit age, but in which this asymmetric Gaussian has 
a dispersion corresponding to the appropriate ltr lower 
or upper value. One interesting difference between the 
histograms and the smoothed distributions can be seen 
when comparing the peaks at ~ 3.5 Gyr in the Z = 0.004 
panel and at ~ 3 Gyr in the Z — 0.008 panel for the 
S99 models. The two peaks are about the same height 
in the histograms, but very different in the smoothed 
version. This difference illustrates how the inclusion of 
uncertainties can alter the significance of peaks. The 
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Fig. 9. — Cluster age histograms derived with the Starburst99 
and GALEV models. The bin size is 10 s yrs. The dotted lines 
indic ate the ages of measured field star bursts (Harris & Zaritskvl 
l200l) . 



Fig. 11. — Smoothed cluster age distributions derived with the 
GALEV models from Monte Carlo simulations (see text for de- 
tails). 
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Fig. 10. — Smoothed cluster age distributions derived with the 
Starburst and GALEV models (see text for details). 



only peaks that appear to reproduce in more than one of 
the smoothed distributions are the ones at ~ 0.5 Gyr, at 
~1.2 Gyr, and at ~2 Gyr. The dotted lines in the Figure 
indicate the position of peaks o bserved in the global sta r 
formation history of the SMC l|Harris fc Zaritskvl F2004). 
Interestingly, the 0.5 and 2 Gyr peaks correspond roughly 
to the peaks observed in star formation history. Are any 
of the peaks seen in the cluster age function statistically 
significant? 

The answer to this question is more complicated than 
simply evaluating the number of clusters of a particular 
age. Our method for estimating ages can result in struc- 
ture in the age distribution even if none is present. If 
the models traverse a large part of color space quickly 
and observational errors scatter clusters within the color 
space, an unrepresentative large fraction of the clusters 



may be close in color to the models that span a small 
age range. To determine how susceptible we are to this 
problem, we simulate our age measurements assuming an 
underlying smoothly exponentially declining cluster age 
function. We make no assumption here as to the cause 
of the exponential decline (whether it reflects a true de- 
cline in cluster formation or the disruption of clusters) . 
We determine colors using both the S99 and GALEV 
models, add uncertainties drawn from our observational 
distribution, and then recover the age distribution using 
our technique. The resulting distributions often show 
some structure, as in the example from the GALEV Z = 
0.004 simulations (Figure ITT|) . To estimate whether the 
structure in these Monte-Carlo simulations is as strong as 
that in the data, we compare the areas under the peaks. 
We automate the measurement of the area and measure 
peak areas for both the data and 1000 simulations (for 
each of the model/metallicity combinations). 

Our principal finding is that no peak is consistently 
significant across the range of models. While for cer- 
tain metallicities, certain peaks are statistically signfi- 
cant (probability of random occurrence < 0.05), without 
knowing exactly which metallicity model is appropriate 
we are unable to conclude that any of the peaks are real. 
In short, both the precision of the colors and the system- 
atic uncertainties inherent in the modeling preclude any 
conclusions regarding the fine structure of the cluster age 
distribution. 

3.2. The Evolution of the Cluster System 

Under the assumption that the cluster formation his- 
tory is the same as that measured for the field stars, a 
comparison of the clust er age distribution and the field 
star formation history ({Harris & Zaritskvl2004D provides 
constraints on the evolution of the cluster system. Any 
variation in the number of clusters normalized by the 
corresponding star formation rate can be interpreted as 
the result of evolution of the cluster system. Various 
factors affect the cluster age distribution, some of which 
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Fig. 12. — The number of clusters of age t per number of stars 
formed over the equivalent time. The two panels show results from 
either Starburst99 (left) or the GALEV (right) models. Within 
each panel the symbols represent the results from either Z = 0.001 
(filled circles), 0.004 (open squares) or 0.008 (crosses). The line in 
the left panel represents a power-law distribution with y oc t~ a ' 1 . 



are observational, such as the loss from the sample of in- 
creasingly fainter or more diffuse clusters with age, and 
some of which are physical, such as the tidal disruption 
or evaporation of clusters. We begin by calculating the 
directly empirical quantity of the number of clusters of 
a given age divided by the number of stars formed at 
that time. We bin our clusters into the same age bins as 
used in the field star formation analysis. This calculation 
provides a star-formation normalized measure of the dif- 
ferential number of clusters as a function of time, which 
approaches dN(t)/dt as the age bin size approaches zero. 
If the population of clusters formed as a fixed fraction of 
the field star formation and if clusters did not evolve, the 
ratio should be constant over time, which is evidently not 
the case (Figure 112)1 . 

The number of clusters currently existing as a func- 
tion of lookback time, dN(t)/dt, is well-approximated 
by the power law function, dN(t)/dt oc t~ , shown in 
the Figure. The best fit slope value for dN(t)/dt as mea- 
sured using both the S99 and GALEV models of different 
metallicities ranges from —1.73 to —2.18, although four 
of the five slopes cluster between values of —2.08 and 
—2.18. In contrast to the details of the cluster age func- 
tion described in §3.1, these results are nearly insensitive 
to the particular model used. To provide some intuition 
for how rapidly the number of clusters in the sample de- 
clines with age we note that for all the best-fit slopes 
more than half the clusters formed are lost less than 10 
Myr after formation, and less than 1% of the clusters 
survive, where survive is simply defined as being part of 
the current sample, past an age of 1 Gyr. 

The observed decline in the number of older clusters 
is due to a combination of selection and physical effects. 
Due to our selection of clusters as resolved stellar sys- 
tems in which we observe stars well down the luminosity 
function, we do not expect strong age-dependent selec- 
tion effects unless the structural properties of the clus- 



ters, such as the central surface density, evolve strongly 
with time. Our conclusion i s similar to tha t reach ed by 
IBoutloukos fc Lamersl {2003) for the lHodgd ltl987l) SMC 
cluster sample. If this claim is correct, then the principal 
cause of the decline in the number of clusters is physical. 

This interpretation is supported by the gen- 
eral agreement we find with previou s studies of 
cluste r disruption both in the SMC ijHunter et alJ 
120031 | B outloukos fc La mersl I20Q3D and other galax- 
ies teo utloTikosfcLarners| I2003IK To be specific, 
IBoutloukos fc Lamersl l)2003D have measured cluster dis- 
ruption in a set of galaxies, including the SMC. In their 
modeling of the problem, they find that they can describe 
the slope of the cluster age function to be (1 — a)/j, 
where a is the slope of the cluster initial mass function 
and 7 describes the power-law mass dependence of the 
disruption time, td{M) oc M 7 . They find that on aver- 
age 7 = 0.57 ± 0.1. From our slope of —2.1 in the SMC, 
we calculate that 7 = 0.48 for an assumed a = 2, which 
is the value of a also assumed by IBoutloukos fc Lamersl 
( 2003). Altern atively, we could adopt the m ean value of 
7 presented bv IBoutloukos fe Lamersl l)2003|) and calcu- 
late a = 2.2, in agreement with the typic al cluster ini- 
tial mass function found in other galaxies l | Zhang fc Fall 
119991: iWhitmore et al.lll999t iBik et al.ll2003|) and within 
the range of 2 to 2 .4 found for the LMC and SMC by 
iHunter et all pOOl . 

One aspect in which we do identify a disagreement 
with previous studies is in the time at which disrup- 
tion do minates over fading fo r the S MC cluster popu- 
lation. IBoutloukos & Lamersl 1)2003)) . using data from 
iHodgd l)1987)K identify a bend in dN(t)/dt, at logt = 8.8, 
that is associated with this transition. We find no un- 
ambiguous flattening of dN(t)/dt at that, or any other, 
time. However, due to uncertainties, we cannot rule out 
a flattening at log (t) < 8.5 for the S99 models, or for 
log(i) < 7 in any of our models. The latter could cor- 
respond to a disruption time scale that agrees with what 
IBoutloukos fc Lamersl i|2003)) find for other galaxies. As 
discussed in §2.1, we have a strong bias against find- 
ing young clusters, and therefore the dN(t)/dt function 
is likely to be even steeper at young ages than what we 
have plotted in FigureEl The lack of a strong flattening 
signature and the uncertainty posed by incompleteness 
preclude us from placing strong constraints on any flat- 
tening at log(i) < 8.5. 

3.3. Investigating the Age Gap 

Th e field star formation history ijHarris fc Zaritskvl 
|2004|) shows two pronounced peaks at about 0.4 and 2.5 
Gyr. As we described above, such peaks are seen in the 
cluster age distribution for certain choices of model type 
and metallicity. While we cannot confirm that these 
peaks exist in the cluster age function, neither can we 
exclude the possibility. Our inability to resolve peaks is 
a consequence of the low precision of our age estimates. 
However, another interesting feature of cluster age func- 
tions, particularly in the LMC, is the presence of an "age 
gap". A similar lull in star format ion is observed in 
the g lobal star formation of the SMC ({Harris fc Zaritskvl 
2004). Can we reach any conclusion about cluster for- 
mation over the timescale of this lull (from 3 to 8 Gyr)? 

In Figure lT3l we show the distribution for cluster ages 
between 1 and 15 Gyr (there are no clusters at measured 
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Fig. 13. — The number of clusters of age t. The two panels 
include the results from either Starburst99 (left) or the GALEV 
models (right) for Z = 0.004. 



Fig. 14. — The number of clusters of age t normalized by the 
expected dissolution function, AiV(i) oc t~ . The two panels 
include the results from either Starburst99 (left) or the GALEV 
models (right) for Z = 0.004. 



ages > 15 Gyr) for both the S99 and GALEV models 
with Z = 0.004. The distributions qualitatively show 
an initial, older epoch of clusters and a recent epoch of 
cluster formation (starting within 3 Gyr of the current 
time), with a relatively quiet time between. The quan- 
titative details, however, are quite different between the 
models. First, the oldest clusters in the S99 models are 
only 8 Gyr old. If this is correct, or if we are incom- 
plete in older clusters, then this pea k would correspon d 
to the ~ 8 Gyr clusters identified bv iRich et al.l {2000). 
On the other hand, if there is a systematic error in the 
entire scale of ages, and the older peak corresponds to 
the oldest LMC clusters which are > 10 Gyr old, then 
perhaps the intermediate peak (currently located at ~ 4 
Gyr) corresponds to the 8 Gyr clusters. Curiously, the 
GALEV models also show an intermediate age peak (at 
~ 6 Gyr), which could correspond to the hypothesized 
8 Gyr population, and these models do identify a truly 
old, > 10 Gyr, population of SMC clusters. In both sets 
of models, the SMC appears to have had an initial clus- 
ter formation epoch followed by relatively little cluster 
formation, with the possible exception of an intermedi- 
ate age peak, punctuated with a recent (< 3 Gyr old) 
episode of cluster formation that is at least as strong as 
the original epoch. 

A different picture of the cluster history is presented 
if we remove the effect of cluster dissolution. In Figure 
1141 we have presented the entire cluster age distribution, 
but we have normalized by the t~ 2A dependence of the 
cluster number. From this Figure it is evident that the 
initial episode of cluster formation dominates the rate of 
cluster formation in the SMC. This Figure rests on the 
assumption that the cluster destruction function remains 
constant in time. From either model, one can conclude 
that the initial flurry of cluster formation (at 11 to 14 
Gyr in the Galev models or 6-8 Gyr in the S99 models) 
did subside, but whether there is an age gap (as seen for 
6 to 11 Gyr in the Galev models) or rather a continuous 
low level of cluster formation (as seen for t < 5 Gyr in the 



S99 models) depends on the model. The more dramatic 
impression of an age gap presented in Figure 1131 is due 
to the fact that the clusters formed in the most recent 
episode of cluster formation have had insufficient time to 
dissolve. 

3.4. Spatial Distribution 

The distribution of clusters of different ages may con- 
strain models of how clusters formed or were destroyed. 
Ivan den Berghl l)1991[) shows that among his sample of 
SMC clusters the young clusters lie along the SMC 
"bar", while the older clusters form more of halo. Al- 
though the_dvnamical reality of the bar has been ques- 
tioned (Zarits kv et al.ll2000f) . the distribution of younger 
and older stars in the SMC follows a similar pattern in 
that the younger stars are preferentially found along the 
north-east/south-west axis while the older stars lie in a 
spheroidal distribution. In Figure [TBI we plot the distri- 
bution of clusters divided into a young (age < 3.5 Gyr) 
and old (age > 3.5 Gyr) samples 1 for S99 and GALEV 
models with Z = 0.004. The characteristics of the distri- 
butions, centroid, ellipticity, and position angle are dif- 
ferent, particularly for the S99 model, which has more 
clusters in the older bin. For example, the right ascension 
centroids of the young and old population arc highly dis- 
crepant (0.942±0.015 and 0.860±0.028 in decimal hours, 
respectively, which corresponds to a difference that is dif- 
ferent than zero at more than 2.6a significance), while 
the declination centroids are consistent (—72.853 ± 0.07 
and— 72.791 ± 0.11 in decimal degrees). 

Because of the poor temporal resolution of our age 
measurements, we are unable to confidently explore other 
age bins, but such a study with higher precision ages may 

1 This div i sion i s quite different than that adopted by 
Ivan den Berghl 119911) . He chose a color cut that for our mod- 
els corresponds to 5.1 X 10 8 and 6.6 X 10* for the S99 and GALEV 
models respectively. Th erefore, his conclusions should not be com- 
pared to our Figure [TBI 
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Fig. 15. — The distribution of young (age < 3.5 Gyr) and old 
(age > 3.5 Gyr) clusters using ages determined with either the S99 
or GALEV models (with Z = 0.004). 



help resolve t he origin of the gas that is hypothesized to 
have infallen ijZaritskv fc Harrisl2004D and the disturbe d 
morphology of the voung stars ijZaritskv et all [2000). 
Correlations among spatial properties and marginally de- 
tected age peaks would also increase the confidence of any 
identified structure in the cluster age function. 

4. CONCLUSIONS 

We use the ima ges from the M agellanic Clouds Pho- 
tometric Survey ijZaritskv et alJ 12002(1 to measure in- 
tegrated U,B,V and / colors for a sample of 195 
stellar clusters with measured structural parameters 
l|Hill fc Zaritskvll2003|) . We estimate cluster ages for 195 
of these using those colors and models of simple stellar 
populations. We conclude that: 

1) Although peaks are visible in the cluster age distri- 
bution, the location of these varies sensitively with the 
choice of metallicity and model type, and similar peaks 
can arise by chance. 

2) We find that the rate of cluster number evolution 
for clusters of age t, after normalizing by the number of 
stars formed at age t, is given by dN(t)/dt oc t~ 2A . 

3) If this decline in the number of clusters with age is 
attributed primarily to physical destruction of the clus- 
ters, we calculate that the disruption timescale at a fixed 
mass is td(M) oc M 0A8 for a cluster initial mass function 
N(M) oc M~ 2 . This result is in agreement with the aver- 
age value measured in a set of nearby galaxies, M°- 57±al 



(|Boutloukos fc Lamersl l2003j) . Alternatively, when we 
adopt the average value for th e disruption timescale from 
iBoutloukos fc Lamersl l|2003D . we calculate that the ini- 
tial mass function N(M) oc M~ 2 - 2 . These values are 
consistent with results from previous studies on other 
galaxies and the SMC. The agreement implies that clus- 
ter dissolution is highly effective and relatively indepen- 
dent of the host galaxy. Only 1% of clusters formed in 
the SMC will survive beyond 1 Gyr. This makes com- 
parison of properties between clusters in the SMC and 
other systems, like the Milky Way, very difficult unless 
ages are known for all the clusters. 

4) There is no evidence of a significant age gap in the 
SMC, although the dominance of the initial epoch of clus- 
ter formation in conjunction with the rapid dissolution of 
clusters produces an observed cluster age function that 
appears to have a formation lull at intermediate ages. 

5) The spatial distribution of the younger and older 
clusters is statistically different, and supports the hy- 
pothesis of a significant accretion or merger event at 
around 3 to 5 Gyr. 

The principal shortcoming of this study is the lack 
of precision in our age measurements that are both the 
result of photometric integrated colors that are suscep- 
tible to large scatter, some of which may be inherent 
to the cluster population, and cluster models that do 
not predict consistent colors for a given age, partly due 
to the lack of metallicity constraints. This problem of 
age and meta l licity degeneracy was also examined by 
Ide Griis et alJ l)2003l) in their study of the systematic un- 
certainties in cluster age determinations in more distant 
galaxies using aperture photometry and stellar popula- 
tion models. The Magellanic Cloud clusters are partic- 
ularly valuable in addressing these issues because ages 
from stellar color-magnitude diagrams can be compared 
to those of integrated colors (§2.3). The current limita- 
tion is the small number of cluster color-magnitude dia- 
grams currently available. To further address the ques- 
tion of systematics, our study can provide guidance as 
to which clusters to target for more precise age measure- 
ments. Such observations benefit both the study of the 
Magellanic Clouds and of more distant systems. 
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